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ABSTRACT 

An overview describing the use of piezoceramic patches in reducing noise in a structural 
acoustics setting is presented. The passive and active contributions due to patches which are 
bonded to an Euler- Bernoulli beam or thin shell are briefly discussed and the results are incor- 
porated into a 2-D structural acoustics model. In this model, an exterior noise source causes 
structural vibrations which in turn lead to interior noise as a result of nonlinear fluid/structure 
coupling mechanisms. Interior sound pressure levels are reduced via patches bonded to the 
flexible boundary (a beam in this case) which generate pure bending moments when an out- 
of-phase voltage is applied. Well-posedness results for the infinite dimensional system are 
discussed and a Galerkin scheme for approximating the system dynamics is outlined. Control 
is implemented by using LQR optimal control theory to calculate gains for the linearized sys- 
tem and then feeding these gains back into the nonlinear system of interest. The effectiveness 
of this strategy for this problem is illustrated in an example. 
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1 Introduction 


The recent development of highly fuel-efficient turboprop and turbofan engines which also 
produce high levels of interior cabin noise (especially at low frequencies) has stimulated a sub- 
stantial effort on the development of a comprehensive active control methodology for interior 
pressure field cavities that have been excited by some primary or external source. In this 
overview paper, we shall discuss recent approaches and preliminary results in the growing ef- 
fort to develop “smart” or “adaptive” material concepts (materials that possess the capability 
for both sensing and actuation are often called “smart” materials) and control strategies for 
such a comprehensive methodology. 

Interior cavity noise in aircraft with turboprop engines is produced primarily through 
(nonlinear) fluid/structure interaction mechanisms. The turboprop blades produce an external 
acoustic pressure field which is converted into mechanical vibrations through fluid/structure 
interactions at the exterior aircraft cabin walls. In turn, these mechanical vibrations produce, 
through interactions of the interior cabin walls with the air in the cabin cavity, pressure waves 
or an interior acoustic pressure field. 

Our discussion here focuses on a time domain state space approach to active or feedback 
control of noise in the interior acoustic cavity. We are especially interested in models and 
methodologies which treat transient dynamics. There is a substantial literature on active 
control of noise in a frequency domain setting (see [18, 21, 24, 26, 28] for some examples and 
further references to both experimental and analytic efforts) as well as a growing literature on 
infinite dimensional state space time domain approaches (e.g. [2, 3, 8, 9, 10]). Earlier efforts 
by most researchers focused on a control methodology implemented through secondary source 
techniques with the input or secondary noise based on feedback of noise levels in the acoustic 
cavity. In this approach, a system of microphones and speakers is strategically placed in the 
interior cavity where one can sense the pressure field (composed of the primary source plus any 
secondary sources present). This information is used as feedback for the actuators or speakers 
which produce a (hopefully) optimally interfering signal (secondary noise) to reduce the total 
noise levels in certain critical zones (related to passenger comfort). Both frequency and time 
domain settings have been used in providing not only “proof of concept” analyses but also in 
designing and and implementing these ideas (to date, mainly in luxury class automobiles). 

More recently, a second approach utilizing smart materials technology has captured the 
attention of investigators. There are a large number of classes of smart materials (e.g. elec- 
trorheological fluids, magnetostrictives, shape memory alloys) but we shall restrict our discus- 
sions in this paper to piezoceramic devices such as piezoceramic patches which, when bonded 
to a structure such as a beam, plate, or curved cylindrical shell, act as an electro-mechanical 
transducer. When excited by an electric field, the patch induces a strain in the material to 
which it is bonded and hence can be employed as an actuator. Moreover, if the host material 
undergoes a deformation (either bending or extension/contraction), this produces a strain in 
the patch which results in a voltage across the patch that is proportional to the strain and 
thereby permits the use of the patch as a mechanical sensor. If constructed and wired with 
proper circuits, these patches can be employed as “self-sensing actuators” [20], thereby provid- 
ing a smart or adaptive material capability for the structure to which the device is bonded or 
in which it is embedded. When combined with a computational adaptive or feedback control 
element, the potential for self-controlled or intelligent structures is enormous. 
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In our presentation and discussions of active control of noise, we shall concentrate on 
actuator aspects of piezoceramics. In the noise suppression example detailed below, we tacitly 
assume that acoustic pressure in the cavity and wall displacements and velocities are sensed 
for feedback. For a complete smart material system, one would use piezoceramic (strain) 
sensors and cavity pressure sensors to construct a state estimator for feedback. 

The motivating example we consider consists of an exterior noise source which is separated 
from an interior cavity by an active wall or plate. This plate transmits noise or vibrations 
from the exterior field to the interior cavity via fluid/structure interactions thus leading to 
the formulation of a system of partial differential equations consisting of an acoustic wave 
equation coupled with elasticity equations for the plate. The control is implemented in the 
example via piezoceramic patches on the plate which are excited in a manner so as to produce 
pure bending moments. It should be noted that the incorporation of the feedback control in 
this manner leads to a system with an unbounded input term (in this case, a system with 
input coefficients involving the Dirac delta “function” and its “derivative”). Experiments are 
being designed and carried out at NASA Langley Research Center in which the interior cavity 
is taken to be cylindrical with a circular active plate to which sectorial piezoceramic patches 
are bonded. 

While the motivating structural acoustics applications are three dimensional in nature, 
many of the theoretical and numerical issues concerning system modeling, the simulation of 
system dynamics, estimation of physical parameters, and the developments of feasible control 
strategies can be studied in 2-D geometries. In this work, we consider a 2-D domain f }(t) 
which is bounded on three sides by hard walls and on the fourth by a flexible beam (see 
Figure 1). A periodic forcing function /, modeling an exterior noise source, causes vibrations 
in the beam which then lead to unwanted interior noise. 



Fig. 1. The 2-D domain . 

This specific problem was chosen since it is a two dimensional slice from a three dimensional 
cylindrical domain which models an experimental apparatus consisting of a rigid cylindrical 
pipe with a clamped aluminum plate at one end. 
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As a 2-D analogue of the plate, the perturbable boundary r 0 (t) (see Figure 1) is modeled 
by a fixed-end Euler-Bernoulli beam having Kelvin- Voigt damping. Bonded to the beam are 
5 pairs of piezoceramic patches which are configured and excited in a manner so as to produce 
pure bending moments (see Figure 2). We reiterate that it is through the excitation of these 
patches that the sound pressure levels are controlled. 

The acoustic response inside the cavity is modeled by a linear wave equation with zero 
normal velocity boundary conditions taken on three walls in order to simulate the rigid walls 
of the experimental pipe. The boundary conditions on the fourth (beam) side of the acous- 
tic cavity result from nonlinear velocity and pressure couplings between the acoustic and 
structural responses (as discussed in [14], these coupling terms are nonlinear since they take 
place along the surface of the vibrating beam). Finally, under the assumption of small beam 
displacements which is inherent in the Euler-Bernoulli theory, the variable domain f l(t) is 
replaced by the fixed domain 0 = [0, a] x [0, as shown in Figure 2. 



Fig. 2. Acoustic cavity with piezoceramic patches creating pure bending moments. 

In terms of the velocity potential <j) (so that p = pj(f)t is the acoustic pressure) and the 
transverse beam displacements w , the strong form of the approximate controlled model for 
the coupled system is then given by 

(j> tt = c 2 A<j> , (x,y) € ft ,t > 0 , 

V<j>-h = 0 , (x,y) € T ,< > 0 , 

V<^(t, x, w(t , x)) • n = w t (t, x ) , 0<x<a,t>0, 

d 2 M 

p Wtt + —- = - pjM t, xM t,x)) + f{t, X) , 0 <x<a, t>0 , (1) 

r\ r\ 

w(t, 0) = ^(t, 0) - w(t, a) = -^(t, a) = 0 , t > 0 , 

(f>(Q,x,y) = (f> 0 (x,y) , w(0,x) = w o {x) 

(f> t (0,x,y) = <l>i(x,y) , u> £ (0,x) = uq(x) 
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(for further details concerning the development of this model, see [14]). Here p, pj and 
c are the beam density, equilibrium density of the atmosphere, and speed of sound in the 
cavity, respectively. The general beam moment M(t,x) consists of an internal component, 
depending on material and geometric properties of the beam and patches, and an external 
component (the control term) which results from the activation of the patches through an 
applied voltage. Specific descriptions of these moments in a variety of settings are given 
in the next section. Finally, the nonlinear coupling between the beam vibrations and the 
interior acoustic field manifests itself in the velocity term V^(t, x, w(t, x)) • n = w t (t,x) and 
the backpressure pf(j} t (t,x,w(t,x)). 


2 Piezoceramic Patch/ Structure Interactions 

As discussed in the last section, control is implemented in the system through the excitation 
of piezoceramic patches which are bonded to the beam. This affects the dynamics of the beam 
in two ways. The first effect is passive and results from the structural changes incurred with 
the bonding of the patches to the structure. In addition to the patch thickness, thefe is a 
nontrivial bonding layer, and both contribute to a moment of inertia which differs from that 
found in regions of the structure not covered with patches. Moreover, the density, Young’s 
modulus and damping coefficient of the glue and patch differ from those of the beam, and 
as a result, these parameters must be modeled as piecewise constants in order to accurately 
match system frequencies (see [17]). The third passive contribution is due to the piezoelectric 
property which dictates that when the patch is subjected to an in-plane strain, a voltage 
proportional to the strain is produced. Hence longitudinal and transverse vibrations in the 
beam lead to the generation of current which provides additional damping in the structure. 
The final (active) contribution from the piezoceramic patches results from the in-plane strains 
which are produced when a voltage is applied. This leads to the generation of external 
moments and forces which enter the equations of motion as external loads. 

The initial part of this section contains a discussion concerning the contributions due to 
patches which are bonded to an Euler-Bernoulli beam. The changes which are necessary for 
extending these arguments to plates and shells are then outlined in the latter part of the 
section with further details given in [16]. 

2.1 Piezoceramic Patch/Beam Interactions 

In the discussion which follows, we consider an Euler-Bernoulli beam of length £, width b 
and thickness h as depicted in Figure 3. The Young’s modulus, mass density (in mass per 
unit volume) and damping coefficient for the homogeneous beam are denoted by Ei>, pi, and 
cd 6, respectively. Bonded to the beam are piezoceramic patches which pan be mounted either 
individually or in pairs as shown in Figures 3 and 4. In the initial discussion concerning 
the contribution due to the patch pairs, it is assumed that both patches have thickness T, 
Young’s modulus E pe , density p pe , and damping coefficient cd pc ■ Moreover, it is assumed that 
the bonding layers for each patch have the same thickness, Young’s modulus, density and 
damping coefficient, and these parameters are denoted by Tu, Ebt, pbe and cpu > respectively. 
We emphasize that these assumptions are made solely for clarity of presentation, and similar 
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results can be obtained in an analogous manner for the more general case in which the patches 
and bonding layers have differing thicknesses and material properties (see, for example, [16]). 

For an Euler-Bernoulli beam having this configuration, force and moment balancing yields 
the strong form of the dynamic equations 


, ,d 2 u dN x A 
P ^ x) dt 2 dx ~ <?r ’ 

d 2 w d 2 M x , 
p ^ dt 2 + dx 2 qn 


drhy 

dx 


(2) 


where N x and M x are the internal force and moment resultants, respectively (see [12, 16]). 
As depicted in Figure 3, w and u denote the transverse and longitudinal displacements, re- 
spectively. The external surface loads q n ,q x and rh y denote normal forces, in-plane forces and 
moments, respectively. For patch pairs with edges at x\ and x 2 , the density of the structure 
is 

p(x) = p h hb + 2b ( puTbt + p pe T) x P e(x) 
where the characteristic function is given by 


{ 1 , X] <x <x 2 

0 , otherwise . 


( 3 ) 



Fig. 3. Cantilever beam with piezoceramic patches. 


A corresponding weak or variational form of the equations can be determined by choos- 
ing V = Hl(T 0 ) x H 2 (T 0 ) for the space of trial functions where r 0 denotes the beam and 
the subscript b again denotes the set of functions which must satisfy the essential boundary 
conditions. Through an energy derivation, one arrives at the variational form 

Jo ^ x ~Bx ~ ~ ° f° r a ll ^i€//f](ro) 

( 4 ) 

Jo + Mx ~dx^ ~ qn< ^ 3 ~ ^ x ~dxJ\ dx = ° * or G ^6 ( r °) 

of the beam equations where N x and M x are external line force and moment resultants. 
As discussed in [16], the surface loads q x and m y of (2) are locally related to the forces and 
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moments N x and M x (which are more natural quantities to use in a weak formulation) through 
the expressions q x = — m y = — ^ t . Global expressions for the specific loads which result 
from the activation of the patches in both the strong and weak formulations are discussed 
later in the section. 



Fig. 4. Strain distribution for the composite structure undergoing bending and extension; (a) patch 
pair } and (b) single patch. 

Internal Moment and Force Resultants 

In order to determine expressions for the internal force and moment resultants N x and 
M xy the patch pair configuration illustrated in Figure 4a is considered first. Because these 
resultants depend upon the stresses and ultimately upon the strains occurring in the structure, 
the description of the resultants begins with a description of the in-plane strains. 

In accordance with the Euler-Bernoulli theory, the strain is assumed to be linear and is 
continuous throughout the combined structure. With e and k denoting the midsurface strain 
and change in curvature, respectively, the strain at an arbitrary point in the beam, bonding 
layer, or patch is given by e = e + kz where z is the distance of the point From the middle 
surface of the beam (see Figure 4a). Because of the differing Young’s moduli and damping 
coefficients in the beam, bonding layer and patch, the stress slopes will differ in the various 
layers. Under the assumption that the stress is proportional to a linear combination of strain 
and strain rate, the stress is given by 

E\>e + cu b e , beam 

a = E b te + cn hl e , bonding layers (5) 

E pe e -(“ co pe e , patches . 

The coefficients co b and CD bt are the Kelvin- Voigt damping coefficients for the beam and 
bonding layer while the coefficient Co pe is taken to be a combination of the Kelvin- Voigt 
damping coefficient for the patch and the damping which results from the production of 
current when the structure vibrates. This latter contribution to the damping results from 
the piezoelectric effect of the patches which dictates that a voltage is produced when the 
patch is subjected to in-plane strains. Under the assumption that the Kelvin-Voigt (material) 
and electrical damping have approximately the same types of effect in the patch, we have 
combined the two into the coefficient co pe which is considered to be unknown and like the other 
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parameters, must ultimately be estimated using data fitting techniques with experimental data 
when considering actual applications. We also point out that the expression (5) can easily 
be generalized to include the possibility of differing material properties in the two patches or 
bonding layers. 

The force and moment resultants are obtained by integrating the stress across the thickness 
of the structure thus yielding the expressions 


II 

f h/ 2 adz 
'—h/2 

II 

r h/2 . 

oz dz 

-h/2 

, regions without patches 

£ 

II 

rh/2 +T bt +T 

/ cr dz 

' — h/2—T b (—T 

II 

o- 

,h/2+T bt +T 
' oz dz 

-h/2-T bi -T 

, regions with patches . 


The substitution of (5) into (6) yields expressions for the resultants in terms of the midsurface 
strain e and change in curvature k. By considering infinitesimal deformations of the middle 
surface, e and k can be related to the longitudinal and transverse displacements u and w 
through the strain-displacement equations 

du d 2 w 

£ = di ’ K = !w 

(see [25], pages 9 and 46). For a beam having two patches bonded to it, the internal (material) 
force and moment resultants are then given by 


N x 


Eh(x)^- + CDh(x) 

ox 


d 2 u 

dxdt 


M x = EI{x) 


d 2 w 


+ C D I(x) 


d 3 w 

dx 2 dt 


( 7 ) 


where 

Eh{x) = E b hb + 26 [E u Tu + E pe T) * pe (x) , c D h(x) = c Db hb + 26 [c Dbt T bl + c Dpe T] Xp e{x) 

El(x) = E b ^ + y [E bl a 3bt + E pe d 3pe ] Xpe (x) , C D I(x) = c D»^y + y [cDlt<*3bl + CD P e a 3pe] XpeC 2 ) • 

Here x P e( x ) again denotes the characteristic function described in (3), and the constants a 3b t 
and a 3pe are given by a 3b t — (h / 2 + T b e) 3 — (h/2) 3 and a 3pe — (h/2 + T b e + T ) 3 — (h/2 + T b e) 3 . 

The substitution of the force and moment resultants in (7) into the dynamic equations 
(2) yields the equations of motion for the combined structure in terms of the transverse and 
longitudinal displacements w and u. As should be expected for a beam containing a pair of 
identical patches which are bonded symmetrically about the middle surface, the differential 
equations (under the first order Euler-Bernoulli assumptions) describing the vibrations in the 
two coordinate directions are uncoupled. 

To see how this differs from the case in which a single patch is bonded to the beam, we 
now consider the case in which a patch of width T is bonded to the beam over the region 
< x < x 2 as shown in Figure 4b. 
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Integrating the stresses through the combined thickness of the structure yields the resultant 
expressions 


, T s du d 2 u d 2 w d 3 w 

N ‘ = Ek ^ai + CDh{z) d^i + E > {x) d? + CD2{x) wdi 

. d 2 w d 3 w du . . d 2 u 

M ' = E, Ww +CDl{x) aMt +E 2 {x) te +CD 2 (x) -Mt- 

The parameters in this case are given by 

Eh{x) = E b hb + b [EhtTu + E pe T) x pc (x) , coh(x) = cr> b hb + 6 [cd&*7m + co pc T] x pe (x) 

^/(x) = T 2 ^ 3 “I" EpcG^pe] ^p e (x) , C/}/(x) = ^ ^ [^/> “ 1 “ pt a 3pe\ Xpe^) 


( 8 ) 


£2(2) — j [^« a 2M + £pe a 2pe] Xpe(*) 


1 c/) 2 (a;) = - [c£) M a 2 w + C0 pe a 2pe ] XpeM 


with ci 3 be and a, 3 pe defined as before and a 2 b( and a 2v e given by a 2 bt = ( h/2 + Tbt ) 2 — (h/2) 2 , 
a 2pe = (h/2 + T be + T) 2 -(h/2 + T bt ) 2 . 

When the force and moment expressions in (8) are substituted into the dynamic equations 
(2), it is apparent that the longitudinal and transverse vibrations are coupled as a result of 
the asymmetry of the structure due to the single patch. This is in contrast to the case when 
patch pairs are bonded to the beam and helps to indicate the, in general, nontrivial effect that 
the patches have on the passive or material properties of the structure. 

External Moment and Force Resultants 

The second contribution from the piezoceramic patches is the generation of external mo- 
ments and forces which results from the converse piezoelectric property that when a voltage 
is applied, in-plane strains are induced in the patch. The magnitude of these induced free 
strains is given by 

e — —V, e — —V> 

c pe 1 J 1 K1 ’ c P e 2 — J' y 2 

where d 3i is a piezoceramic strain constant, and V\ and V 2 are the voltages into the two 
patches in the pair. We point out that when a voltage is applied to a free patch with edge 
coordinates x\ and x 2 , the point x = (a?i + x 2 )/2 will not move whereas the symmetric points 
on either side wifi move an equal amount in opposite directions. This motivates the use of 
the indicator function in several of the following definitions. 

The stresses due to the excitation of the patches are given by 


(&x)pe 1 — E pe e pe , , (®x)pe 2 — E pe € 


pe 2 


with the negative signs resulting from conservation of forces when balancing the material and 
induced stresses in the patch. 

The integration of these stresses through the thickness of the patches yields the expressions 


(M x ) pe = [(M x ) 

pe 1 + (M X ) 

pe 2 ] Xpe(z) 

(JV.)„ = [(JVJ P e l + (A*) Pe 2 ] Xpe(^) 5 ’l, 2 (®) 


( 9 ) 


M 

I 
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where 


(M x ) pei = - l -E ve bdv{h + 2 T u + T)Vi 
(•^x)pe i — Ep e bd>3 \V\ 


(Mc)pe 2 = - E pe bd 3 \{h + 2 The + T )\ 2 
{N x )pe 2 = — E pe bd 3 \V 2 


for the external moments and forces generated by the activation of the patches. The presence 
of the indicator function 


SxM 


' i 

< 0 
-1 


5 


5 




x < (xj + x 2 )/2 
x = (rci + x 2 )/2 
x > (xi + x 2 )/2 


( 10 ) 


results from the fact that for homogeneous patches having uniform thickness, opposite but 
equal strains are generated about the point x = {x\ + x 2 )/2. 

These expressions can then be substituted directly into the weak equations (4) as loads on 
the beam (with q n = 0 and N x = ( N x ) pe , M x — ( M x ) pe ). In order to determine the patch loads 
for the strong form of the beam equations, the corresponding surface moments and forces are 
found via the relationships 


q x = -Si, i{x) 


djNxjpe 

dx 


rh y 


djAQpe 

dx 


and these latter values are used in (2). We point out that this results in the need to differentiate 
across discontinuities in characteristic and indicator functions (once for the force and twice 
for the moment) whereas this problem is avoided in the weak formulation since the derivatives 
are transferred on the test functions. In fact, the effect of the characteristic functions in the 
latter case is to simply restrict the integrals to the region covered by the patches. 

The general moments in the beam component of the structural acoustic system (1) can 
now be described in terms of the internal and external moments just discussed. By combining 
both the passive and active contributions due to a single pair of patches which are excited 
out-of-phase, the general moment is given by 


M = M x - (M x ) pe 


where the internal and external moments are 

d 2 w . . d 3 w 

M, = El(x) w + cdI(x) T x Wt 

( M x ) pe = E pe bd 3 i(h + 2Tu + T)Vxpe(x) = )C B Vxpe(x) 


as given in (7) and (9), respectively (the latter expression is obtained by taking V — V\ = —V 2 
in (9)). We emphasize that the out-of-phase excitation of the patches produces pure bending 
moments and hence only transverse vibrations are present in the beam response. 

For a system in which s pairs of patches are bonded to a beam and are excited out-of-phase, 
the beam component of the system (1) has the form 


d 2 w d 2 


d 2 w 


d 3 w \ 


d 2 


p{x) W + w [ EI{x) w + CdI(x) wNJ + P/Mw) = f + E«W)^x w (*) 
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where Xpe,( x ) denotes the characteristic function over the i th patch pair and Ui(t ) is the voltage 
into the i th pair. The parameters El and cpI are given by 



N h 3 b ’ 2b r , 1 / X 

CDi(x) — C£> [, + / ^Db(, ^3M, d" c Z?pe, a 3pe,J Xpe.t^j 


1=1 


while the patch parameters are given by Kf = E^Rdzi^h + 2 Tbt, + 7)) (in these definitions, 
the bonding layers and patches in the i lh pair are considered to have thickness Tu, and 
Ti, respectively). We note that the discontinuous parameters p,EI,C£>I and Kf lead to 
second derivatives of characteristic functions which causes difficulties in the strong form of 
the equations. The transfer of these derivatives onto test functions eliminates these problems 
in the weak form of the equations and is one motivation for using the weak form of the system 
equations as discussed in the next section. 


2.2 Patch Contributions to Plate and Shell Dynamics 

In the first part of this section, the contributions from piezoceramic patches to the longitudinal 
and transverse vibrations of an Euler-Bernoulli beam were examined. It was noted that the 
patch contributions could be categorized into two types; the first resulted from the structural 
changes incurred when the patches were bonded to the beam while the second effect was due 
to the activation of the patches when a voltage was applied. These same types of effects 
result when piezoceramic patches are bonded to more complex structures such as thin plates 
or shells. 

The motion of a plate differs from that of a beam in that two sets of longitudinal motion 
are present with the stretching in one coordinate direction related to the contraction in the 
other through the Poisson ratio v. In thin shells, the transverse and longitudinal vibrations 
are coupled due to the underlying curvature of the structure. However, once the underlying 
dynamic equations in terms of the force and moment resultants are known, the effects due 
to the presence and activation of the piezoceramic patches can be determined in a manner 
analogous to that discussed above for thin beams (see [16]). 

To illustrate, we consider a thin circular cylindrical shell of radius R , thickness h and 
having the axial coordinate x as shown in Figure 5. As in the beam discussion, the variable z 
measures the distance of a point on the structure from the corresponding point on the middle 
surface ( z = 0) along the normal to the middle surface. 
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Fig. 5. The thin cylindrical shell. 

As discussed in [16], the infinitesimal strain relationships for a cylindrical shell are 


e x = e x + zk x 


1 


eg = 


1 + z/ R 


(e e + zk 9 ) , 7 x e = 


1 


l + z/R 


+ z (' + Jr) ti 


where e x and eg are normal strains at an arbitrary point within the cylindrical shell and ^ x g is 
the shear strain. Here e x , eg and e x g are the normal and shear strains in the middle surface and 
k x , K$ and t are the midsurface changes in curvature and midsurface twist (see [23], page 8). 

In terms of the axial, tangential and radial displacements u, w and w, respectively, the 
expressions for the midsurface strains and changes in curvature for the cylindrical shell are 




Kx = 


du 

dx 

d 2 w 


1 dv w dv 1 du 

£e= Rd9 + R ’ £x6= dx + R'de 

1 d 2 w 1 dv _ 2 d 2 w 2 dv 

Ke = + R^de ’ T = ~Rdxd6 + Rd^ ' 


If a generalized Hooke’s law in which stress is assumed to be proportional to a linear combi- 
nation of strain and strain rate is used as the constitutive relation, the stresses in the shell 
are given by 

<J X = - + v » e6 ) + i CP ' *" 2 (g* + v »te) 

l - v; i - 

R3 / , \ . ^Ds / • . • \ 

ere = Aeg + v,e x ) + -{eg + v„e x ) 

1 - v 2 I -v 2 . 


Es . CDs 

CTx6 — &6X — n( 3 , \7 X$ T 0 / ! , 

2(1 + v s ) 2(1 + Vs) 

where <j x and ag are normal stresses and a x g and a g x are tangential shear stresses. The 
constants E s , v s and cd s are the Young’s modulus, Poisson ratio, and damping coefficient for 
the shell. Similar relations are found in the bonding layers and patches (see (5) for analogous 
expressions for the beam). 
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The internal or material moment and force resultants are obtained by integrating the 
stresses across the thickness of the structure. For patches having thickness T and bonding 
layers of thickness Tm, this yields the expressions 


N x 

N x9 


M x 

M. 


x$ 


rh/ 2 +T bt +l 

J-h/2-T be -: 

J—h/2—Tbt— 


G X 

(\ + ^]dz 

' N e 

rh/2+T bt +T 

<r$ 

&x6 

V RJ 

Ng x 

J-h/2-T bt -T 

&6x 


(T x 

( 1 + t;! z dz , 

' Me 

rh/2+T bt +T 

ere 


V RJ 

Mq x 

J-h/2-T bt -T 

ve x 


dz 


zdz 


(ii) 


in regions of the structure covered by the patches with similar expressions in those region of 
the structure consisting solely of shell material (the limits of integration in this latter case are 
—hf 2 and hj 2). Explicit descriptions for these internal moment and force resultants can be 
found in [16]. 

In a shell which is excited by the activation of piezoceramic patches, the external moments 
and forces are due to the in-plane strains 


°V e \ ~ ( e x)pe i — ( e fl)pei — ' rjn K 


Z V e 2 ~ ( e z)pe 2 — ( e e)pe 2 — rp V? 


which result from the input of the voltages V\ and Vi into the outer and inner patches. The 
resulting external stresses are given by 


(Vx)pe\ — { cr e)pe i — 


E i 


1 - z/, 


-pe i 


(° r r)pe 2 — ( a o)pc - 2 — 


E o 


1 — 1/2 


-"P e 2 


For a patch with bounding values X\ , X 2 , 0\ and 0 2 the total external line moments and 
forces are 

(M x )p e = [(A/ I ) p e 1 + ( A/ X )p e2 ] X pe (x, 0) 


(M e ) pe = {(M e ) pe x + (M e ) 

V e 2 ] Xpe( x i 0) 

(N x )pe = [(A(r)pei + (^x)pe 2 ] Xpe( x , 0)Si, 2 (x)Si i2 (O) 

(iV^)pe = [(A^pe, + (Ne) P e 2 ) Xpe( x ,O)Sl,2(x)S it2 (0) 


( 12 ) 


where the indicator function S'j^x) is defined in (10) (with a similar definition for S\ >2 (0)) 
and 

, m f 1 , X! < x < x 2 , 0, < 0 < 0 2 

Xpt{ x iO) — ( o , otherwise. 

The individual patch moments are obtained by integrating the external stress distribution 
through the thickness of the patches in the same manner used in the beam analysis (see [16]). 
In order to obtain a strong form of the equations of motion, force and moment balancing 
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im m 



can be used to obtain Donnell-Mushtari shell equations 


Rp( ' e ' dt 1 R dx 


dN„ = _ r 3(^ i±Si ^ x )s hl (6) 


Rp(x, 9) 
Rp{x,9)^ - R? 2M * 


&v dN e 
dt 2 dO 


R 


89 

N*1 

dx 


d(N e ) pe 

89 


s hl (x)ih, 2 («) 


dt 1 


dx 2 


1 d*M, „3 2 iW rf , „ D . „9 2 (M t ) p , 
TW- - 2 ‘ teW + N> = ~ R ~ dT~ 


1 <9(A/fl) pe 
R 89 


(see [16, 23] for a more detailed derivation of these equations as well as a discussion concerning 
the assumptions that are made in obtaining this and other forms of the equations of motion for 
a thin shell). The contributions due to the patches are incorporated in the internal moments 
and forces (11), the external moments and forces (12), and the variable density p(x,9). 


3 Weak Form and Well-Posedness of the Structural 
Acoustics Model 

As discussed in the last two sections, the incorporation of the piezoceramic patch contributions 
into the strong form of the modeling system equations leads to first and second derivatives of 
characteristic functions since both the internal and external moments contain discontinuities 
at the edges of the patches. This yields an unbounded control input operator and leads to 
difficulties when approximating the dynamics of the coupled system. To avoid these difficulties, 
it is advantageous to formulate the problem in weak or variational form (the use of the 
variational form also permits the use of basis functions having less smoothness than those 
used when approximating the solution to the strong form of the equations). 


3.1 Weak Form of the System Equations 

The state for the second-order form of the 2-D structural acoustics problem is taken to be 
z = in the Hilbert space H = L 2 (fi) x L 2 (T 0 ) with the energy inner product 

The choice of the space Z 2 (fl), defined as the quotient of L 2 (fl) over the constant functions, 
results from the fact that the potentials are determined only up to a constant. 

To provide a class of functions which are considered when defining a variational form 
of the problem, we also define the Hilbert space V = H X (H) x H${To) where H\Cl) is 
the quotient space of H l over the constant functions and Hq(To) is given by /^(Fo) = 
{ 0 g H 2 (T 0 ) : i>(x) — 0'(x) = 0 at x = 0, a}. The V inner product is taken as (here and below 
we use the notation D = ^) 

pjV<f> • V(du> + f EID 2 wD 2 rjd'y . 
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As discussed in [14], integration in combination with the use of Green’s theorem then yields 
the nonlinear first-order variational form 

/ + 

Ja c z 

+ 

+ 


for all (£, 7]) in V (here Xpe,(%) denotes the characteristic function over the i tfx patch). We note 
that the nonlinear coupling term can be written as <f> t (t , x, w(t, x)) = <f> t (t , x, 0) + ^6f (i , x, x)) 
where <^(£, x, w(t, x)) = <^ t (i,x,itf(f,x)) — <^(£,x,0). We will make use of this decomposition 
in the abstract formulation of the nonlinear system as a perturbation of a linearized system 
in our discussion below. Again, a more complete discussion and motivation concerning the 
formulation of the first-order system in weak form is given in [14], 

We point out that in this variational form the derivatives have been transferred from the 
plate and patch moments onto the test functions. This eliminates the problem of having to 
approximate the derivatives of the characteristic function and the Dirac delta as is the case 
with the strong form of the equations. 

The system (13) can be formally approximated by replacing the state variables by their 
finite dimensional approximations and constructing the resulting matrix system. Hence it is 
in a form which is suitable for use in applications. In order to discuss the well-posedness of 
the model, however, it is advantageous to pose the problem in terms of sesquilinear forms and 
the bounded operators which they define, and this is the subject of the rest of the section. 


/ pw ti rjd 7 

^r 0 

f prV(j> • + f EID 2 wD 2 rjd 7 

J 0 Jr Q 

f {c D ID 2 w t D 2 T) + p f [<f>t(w)ii - Wtfljd'y 
‘'To 

/ i2 K f u ^ t )xv^x x ) D2r i d i + L fv d 7 

,- = | ''To 


(13) 


3.2 Abstract First-Order Formulation 

As motivated by the theoretical results in [3, 4, 6, 15], we consider (see [29] for basic definitions 
and fundamental functional analysis theory) the Gelfand triple V *— ► H ~ H * <— ► V* with pivot 
space H and define sesquilinear forms eq : V x V — ► 1R, i = 1,2 by 


<ti(4>, 4*) = / pjV<f> ■ V£du> + / EID 2 wD 2 T]d'y , 


< t 2 ( 4 *, ) = / {c d ID 2 wD 2 t] + pj(<f>t] - w()}dy 

Jr 0 ....... 

where $ = (<£, w) and ^ = (£, 77 ) are in V. 

As detailed for a similar problem in [15], it is straightforward to show that with these 
definitions, <7i and <r 2 are bounded (there exist c\ and c 2 such that |<7i(<J>, ty)] < ci |4>|v|'I'|i/ 
and |<x 2 (<I>, $)| < c 2 |4>|v|^|v), &i is V-elliptic and cr 2 is H-semielliptic (there exist c > 0 and 
b > 0 such that Re <7i(<I>, $) > c|4>|v, and Re<7 2 ($,$) > for all 4> € V) and that a\ is 
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symmetric (< 7 i($, $) = $) for all $, € V). As a result of the boundedness, we can 

define operators A\, A 2 € £(V, V*) by 

for i = 1,2. 

To account for the control contributions, we let U denote the Hilbert space containing 
the control inputs (U = IR in our structural acoustics example), and we define the control 
operator B € £(U,V*) by 

(Bu,V) v . y = [ ^JCfuiXpe.WD 2 ^ 

for ^ e V, where (v)v*,v ' s ^ ie usua l duality pairing. Finally, letting F = ( 0,f/pb ) and 

G(z,z t ) = (0 ,-PfMw)) where again, } t (w) = j> t (t,x,w{t,x)) = &(«,*, w(*i*)) - 0) 

denotes the nonlinear perturbation to the linear coupling term, we can write the control system 
in weak or variational form 

(z ti (t), 'k)v i y + 02(~t(*)> v k) + < 7 j(z(*), '!') = { Bu(t ) + F(t) + G(z(t),z t (t)), ( 14 ) 

for <]> in V. This then yields the system 

z tt {t) + A 2 z t {t) + A^z(t) = Bu(t ) + F(t) + G(z(t),z t (t)) 


in V*. 

To apply infinite dimensional control results for periodic forcing functions to this problem, 
it is advantageous to write the system in first-order form. This is accomplished by defining, 
the product spaces 7 H = V x H and V = V x V with the norms 

k$> n 2 n = i^iv + i^i// ) 

\(*,m = m 2 v + m 2 v 

We point out that V •— ♦ H ~ H* *-> V* again forms a Gelfand triple. 

The sesquilinear form a : V x V — ► IR is then defined by 

cr(0, x ) = <r((T, A), (*, ¥)) = - (A, *) v + a,( T, ¥) + cr 2 (A, ¥) 
where x = ($> ’P) and 0 = (T, A). 

For the state Z(t) = (z(t),z t (t)) in 7i, we can subsequently write the system in the first- 
order variational form 

(Z t (t), x) v *,v + *(2(0, X) = {Bu(t) + T{t) + Q{Z{t)\ x) v *,v ( 15 ) 

where T{t) = (0 ,F(t)), G(Z(t)) = (0, G(z(<), z t {t))) and Bu(t) = (0 ,Bu(t)). As usual, the 
relation (15) must hold for all x € V. Finally, the weak form (15) is formally equivalent to 
the system 

Z f (t) = AZ(t) + C{t y Z(t)) (16) 
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in V* where 


(17) 


and 


C{t, Z{t)) = Bu(t) + Fit) + G{Z(t)) 

dom A = {0 = (T, A) e H : A € V, AjT + A a A € H} 

[0 I ] 

A = 

-A x -A 2 


(18) 


3.3 Model Well-Posedness 

In the previous discussion, the weak form of the coupled structural acoustic equations was 
written as an abstract first-order semilinear initial value problem with a state in / H. The 
nonlinear forcing term C(t,Z(t)) = Bu(t) -f- !F(t) + however lies in V* rather than 

H since the control term B € C(U, V*) defines the product space control term Bu(t) = 
(0, Bu(t)) € {0} x V* C V x V* = V*. Hence the standard theory for abstract semilinear 
Cauchy problems does not apply directly, and the first step in the following discussion is the 
outline for arguments which can be used to extend the operator A to a space where the theory 
does apply. A more extensive discussion concerning the well-posedness of a linear problem of 
this type can be found in [15] and details for the following arguments can be found in that 
work. 

The first step in determining the well-posedness of the system model is to argue that 
A generates a Co-semigroup on Ti. As noted earlier, the sesquilinfiar form <J\ is V-elliptic, 
continuous and symmetric while a 2 is continuous and H-semielliptic. From the Lumer-Philips 
theorem (with further arguments found in [1] and pages 82-84 of [4]) this then implies that 
the operator A defined in (18) generates a Co-semigroup on the state space Ti. Moreover, the 
semigroup satisfies the exponential bound |T(f)| < e ut for t > 0 (where in fact, u = 0 due to 
the fact that A is dissipative as shown in [4]). 

Since Bu(t ) lies in V* rather than H, the next step is to extend the semigroup T (t ) on H 
to a semigroup T(t) on a larger space W* D {0} x V* so as to be compatible with the forcing 
term (this is accomplished using “extrapolation space” ideas and arguments similar to those 
presented in [6, 7, 22]). 

As detailed in [15], the space of interest is defined in terms of dom *4* where 

dom A* = {x = ($, * ) € n\ * e V , - A* 2 V € H) 


A* X = 


a;$ 


-* ) . 

-a;* ) 


Specifically, the space W = [dom A *] is taken to be dom A* with the inner product 


for some arbitrary but fixed A 0 with A 0 > w (recall that the original solution semigroup 
satisfies the bound \T(t)\ < e wt ). As proven in [7], the resulting W norm is equivalent to 
the graph norm corresponding to A*. Moreover, we have that {0} x V* C W* = [dom «4*]* 
(see [15] for details). 
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From the definition of A* and the equivalence of the W norm with the graph norm corre- 
sponding to A *, we can define AQ £ W* by 

(.46) (x) = (&,A'x) n 

for all 0 £ H, X € W. With this definition and the Riesz representation theorem, it is shown 
in [15] that A is an extension of the original operator A from dom.4 C H to all of Ti. Finally, 
as proven in [7], the operator A is the infinitesimal generator of a Co-semigroup T(t) on W* 
which is an extension of T{t) from H to VV*. 

In the corresponding linear problem, under reasonable regularity conditions on t u(t ) 

and t »-» F(f), one can immediately argue the existence of a unique strong solution to the 
system in terms of the extended semigroup 7~(/). For the semilinear problem of interest, 
however, the nonlinear nonhomogeneous terms must satisfy certain continuity criteria in order 
to obtain similar results. For example, if we let X denote the reflexive Banach space W* and 
assume that C : [0,T] x X — ► X defined in (17) is continuous in t on [0, T] and uniformly 
Lipschitz continuous on X , then the integral equation 

z (t) = tm + /; nt - .) ( Bu(s) + F( ° + G(2(s)) ) * ( 19 ) 

is well-defined for Bu + F + G(Z) G L 2 ((0, T), V*). Moreover, for Z(0) = Z 0 , the solution Z(t) 
of (19) is a unique mild solution to (16) (see Theorem 1.2, page 184 of [27]). In addition, if 
C : [0, T] x X — » X is Lipschitz continuous in both variables, then it follows from Theorem 1.6, 
page 189 of [27] that (19) provides the strong solution to (16) interpreted in the W* sense. 

The required continuity of the nonhomogeneous terms Bu and F is demonstrated in [15] 
and hence the remaining question concerns the Lipschitz continuity of the nonlinear coupling 
term G(z,z t ) = (0, — P/<f>t{w))- If we assume that the input terms F and Bu are sufficiently 
smooth so as to assure the necessary continuity in G(z,z t ), then our open loop nonlinear 
system is well-posed. 

3*4 Well-Posedness of the Closed Loop System 

The arguments leading to the well-posedness results for the linear and nonlinear open loop 
models can also be extended to the closed loop systems which result when the gains determined 
for a corresponding LQR problem are fed back into the system. In determining these gains, 
the perturbing force T is assumed to be periodic (this is a reasonable assumption since JF 
models the exterior noise which in this problem is generated by the revolution of turboprop 
or turbofan blades). 

Discussing first the linearized problem, the periodic LQR problem consists of finding u G 
L 2 (0,r;[7) which minimizes a quadratic cost functional of the form 

J(u)= l - [ T {(QZ(t),Z(t))n + {Ru(t),u(t))u} dt 
2^ J o 

subject to Z t (t ) = AZ(t) + Bu(t) + F(t) with Z( 0) = Z(t). Since Z = (</>, w, <f> t , w t ), the 
operator Q can be chosen so as to emphasize the minimization of particular state variables as 
well as to create windows that can be used to decrease state variations of certain frequencies. 
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The control space U is taken to be IR* if s patches are used in the model, and it is assumed 
that the operator R € C,{U) is an s x s diagonal matrix where r„ > 0,i = is the 

weight on the controlling voltage into the i th patch. In the case that B is bounded on ? i, a 
complete feedback theory for this periodic problem can be given as discussed in [19]. This 
theory can be extended to also include the case of unbounded B, i.e., B € £((/, V*), of interest 
here (see [5]). Under usual stabilizability and detectability assumptions on the system as well 
as standard assumptions on (2, the optimal control is given by 

u(t) = -R- l B*[UZ(t) - r(t)] (20) 

where II € £(V*,V) is the unique nonnegative self-adjoint solution of the algebraic Riccati 
equation 

A*n + YlA-UBR-'B*n + Q = 0 . (21) 

Here r is the unique r-periodic solution of 

r(t) = -( A * - UBR-'B*)r(t) + ILF(f) (22) 

and the optimal trajectory Z is the solution of 

Z(t) = (A- BR~'B*U)Z(t) + BR~ l B*r{t) + F{t) . (23) 

As discussed in [5] for the case when B € C(U, V*), one also finds that the operator 
A— BR~^B*Xl generates an exponentially stable Co-semigroup «!>(£) on the state space Ti. From 
Corollary 10.6, page 41 of [27], this implies that A* — UBR~ 1 B* generates the corresponding 
adjoint semigroup S*(t) on ~ fi. The semigroup S(t) can then be extended through the 
extrapolation space techniques just discussed to a larger space VV* 3 {0} x V *, and with 
reasonable regularity assumptions on f 1-4 F'(t), this implies the existence of solutions to the 
tracking equation (22) and closed loop system (23) for r(0) = r 0 and Z( 0) = Z 0 . 

As discussed in greater detail in the next section where the corresponding finite dimensional 
control problem is considered, an effective strategy for controlling the original nonlinear system 
is to determine the gains for the linearized model and feed these back into the nonlinear system. 
This then yields the nonlinear closed loop system 

Z{t ) = (A - BR- x B*W)Z{t) + BR~'B*r(t) + T(t) + Q{{Z{t)) 

where again, A — BR~'B* II generates the Co-semigroup S(t ) which can then be extended 
to VV*. With the assumption that the input term F is sufficiently smooth so as to assure 
the necessary continuity in nonhomogeneous terms, the closed loop nonlinear system is also 
well-posed. 
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4 System Approximation and the Finite Dimensional 
Control Problem 

The discussion thus far has centered around the infinite dimensional model for the structural 
acoustic system as well as issues concerning its well-posedness. However, in order to develop 
viable schemes for approximating the nonlinear system dynamics, estimating physical param- 
eters, and determining control gains, appropriate finite dimensional approximations to the 
state variables w and <f> must be developed. For reasons discussed in [3], a Galerkin scheme 
was chosen and the potential and beam displacement were discretized in terms of spline and 
spectral expansions, respectively. 


4.1 System Approximation 

A tensored Legendre basis was used for the discretization of the acoustic velocity potential. 
Letting P t a (x) and Pf(y) denote the standard Legendre polynomials that have been scaled by 
transformation to the intervals [0, a] and [0,4], respectively, the basis functions {B™} for the 
cavity were then defined as 

B%{x,y) = P?{x)P*(y) for i = 0, 1, • • • ,m„ j = 0, 1, • • • ,m y , i + j ^ 0 , 

where m = ( m x + 1) • (m y + 1) — 1. The condition i + j ^ 0 eliminates the constant function 
thus guaranteeing that the set of functions is suitable as a basis for the quotient space. The 
m dimensional cavity approximating subspace is taken to be H™ = span and the 

approximate cavity solution is given by 

771 

1 = 1 

771 y m x 

= EE %m p ?( x ) p ‘(y) ■ 

j =0 »=o 

o 

Cubic splines were used as a basis for the beam displacement since they satisfy the smooth- 
ness requirement as well as being easily implemented when adapting to the fixed-end boundary 
conditions and patch discretizations. Letting {B ™ denote the cubic splines which have 
been modified to satisfy the boundary conditions (see [3, 14] for details), the corresponding 
n — 1 dimensional beam approximating subspace is given by H£ = span{B t n }^ =1 and the 
approximate beam solution is taken to be 

i=l 

The approximating state space was then taken to be H N = H™ x Hg where N = m + n— 1, 
and the product space for the first order system is / H N = H N x H N . By restricting the infinite 
dimensional system (13) to 7i N x H N , one obtains the nonlinear finite dimensional system 

M N y N (t ) = A N ( 2 y N (t )) + B N u(t) + F N (t) 

M N y N ( 0) = y» 
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or equivalently 


y N (t) = A N (i y N (t )) + B N u(t) + F N (t) 

y N ( o) = y? • 


(24) 


Explicit descriptions of the mass and stiffness operators M N and A N {y N {t)) as well as detailed 
definitions of the control matrix B N and the force vector F N (t) can be found in [3, 13]. The 
vector y N (t) = ■ ■ • , w?(t), • • • , «£_,(<), ‘ * * » ^m(0> ^(0* ' ’ ' » wJL,(«)) r 

contains the 2N x 1 approximate state coefficients while u(t) = («i(Q, • • • u s (t)) T contains the 
s control variables. As detailed in [14], the nonlinearity in the operator A N (y N (t)^ manifests 
itself in the dependence of the operator on the unknown coefficients W*)}- 


4.2 The Finite Dimensional Control Problem 

Due to the nonlinearity in the infinite dimensional system (13) and hence the finite dimen- 
sional matrix system (24), LQR feedback control results for problems with periodic forcing 
terms can not be directly applied as there were in [3]. Instead, the following strategy was 
adopted. The infinite dimensional system was linearized by replacing the nonlinear cou- 
pling term x, w(t, x)) by its linear component <f> t (t,x, 0) (this is equivalent to taking 
G(z(t), z t (t)) = 0 in (14) or ( Q{Z{t )) = 0 in (15) or (16)). This linearization is motivated by 
the assumption of small beam displacements which is inherent in the Euler-Bernoulli theory 
(for physically reasonable input forces, the beam displacements are of the order I0 _5 m for 
the geometries of interest). The feedback gains for this approximate linearized system were 
calculated from a periodic LQR theory (see [3]) and were then fed back into the nonlinear 
problem to create a stable nonlinear closed loop control system. 

To illustrate this control strategy, the LQR theory for problems with periodic input terms 
is briefly outlined. The resulting gains are then applied to the nonlinear problem of interest 
with the results being illustrated in an example. 

Linear Periodic Control Problem 


As discussed in [3], the approximation of the nonlinear coupling term <f> t (t,x,w(t,x)) by 
its linear component, and the projection of the resulting system into the finite dimensional 
subspace x 7i N yields the linear finite dimensional Cauchy equation 


y N {t) = A N y N (t ) + B N u(t) + F N {t ) 
»"(«) = Vo 


(25) 


(this system can also be obtained by restricting the infinite dimensional system (14) with 
G(z(t), z t {t)) = 0 to x H n ). The components of the linear stiffness matrix can be found 

in [3]- 

The periodic finite dimensional control problem is then to find u € L 2 ( 0, r) which minimizes 


J N ( «) = \£ {< «V(t ). + <*«(«). «(())«•} dt , N = m + n - 1 
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where y ^ solves (25), t is the period, R is an s x s diagonal matrix and ra > 0, i 1,‘ ■ ■ ,s 
is the weight or penalty on the controlling voltage into the i th patch. 

The nonnegative definite matrix Q N is chosen in a manner so as to to emphasize the 
minimization of particular state variables. From energy considerations as discussed in [3], an 
appropriate choice for Q N in this case is 

Q n = M n V 

where M N is the mass matrix, and the diagonal matrix V is given by 

V = diag • 

Here I k ,k = m,n — 1 , denotes a k x k identity and the parameters d{ are chosen to enhance 
stability and performance of the feedback. 

The optimal control is then given by 

u N (t) = R~ 1 (B n ) t [r*(t) - n N y N (t)] 


where is the solution to the algebraic Riccati equation 


(A N ) T n N + u n a n - n n b n r~' (b n ) t u n + q n = o . (26) 


For the regulator problem with periodic forcing function F N (t), r^(f) must satisfy the linear 
differential equation 



r"(0) = r w (r) 


B r ''R-'(B !l ) r n r -'\ T r"(t) + n N F N {t) 


(27) 


while the optimal trajectory is the solution to the linear differential equation 

y»(t) = [A N - B N R-\B N ) T n N ] y N (t ) + B N R-\B N ) T r N (t) + F N (t) 
y N { 0) = v n {t) ■ 

The finite dimensional optimal control, Riccati solution, tracking equation and closed loop 
system can be compared with the original infinite dimensional relations given in (20), (21), 
(22) and (23), respectively. In order to guarantee the convergence — * II, r N —* r, and hence 

the convergence of u N — * u, it is sufficient to impose various conditions on the original and ap- 
proximate systems. These hypotheses include convergence requirements for the uncontrolled 
problem as well as the requirement that the approximation systems preserve stabilizability 
and detectability margins uniformly. A fully developed theory (see [5]) is available for the 
case when T = 0 (in this case the tracking variable r does not appear in the solution) even 
when B is unbounded. However, the theory in [5] requires strong damping in the second-order 
system whereas the only damping in our system is the strong Kelvin- Voigt damping in the 
beam (damping in the cavity was omitted due to the relatively small dimensions involved). 
Although the convergence theory of [5] does not directly apply here, numerical tests indi- 
cate that convergence is obtained even though this system contains only weak or boundary 
damping. 
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Nonlinear Control Problem 

To extend these results to the nonlinear system of interest, the linear gains were calculated 
and fed back into the nonlinear system (24), thus yielding the suboptimal control 

u N {t ) = R-\B n ) t [r w (t) - nV(<)] 

and the closed loop system 

y N {t) = A N (y N {t)) - B N R-'{B N ) T n N y N (t) + B N R~'(B N ) T r N (t) + F N (t) 

V N { 0 ) = y N (r) ■ 

The Riccati matrix fl^ and tracking vector r N (t) are solutions to (26) and (27) which arise 
when formulating the corresponding LQR problem. 

Example: Nonlinear Control 

To illustrate the dynamics and effects of feedback control on a nonlinear system mod- 
eling a 2-D analogue of a 3-D experimental setup, a .6 m by 1 m cavity with a flexible 
beam at one end was considered (see Figure 6). The beam was assumed to have width 
and thickness .1 m and .005 m, respectively, and the Young’s modulus and beam density were 
taken to be E = 7.1 x 10'° N/m 2 and pb = 2700 kg/m 3 . This yielded the stiffness parameter 
El = 73.96 Nm 2 and linear mass density p = 1.35 kg/m. The damping parameter for the 
beam was chosen to be c^I = .001 kg m 3 / sec. The speed of sound and atmospheric density 
inside the cavity were taken to be c = 343 m/sec and pj — 1.21 kg/m 3 , respectively. 

1 


0 .6 

FlG. 6. Acoustic cavity with one centered 1/2 length patch. 

Several forcing functions modeling uniform (in space) periodic exterior sound sources were 
considered. In this example, the forcing function was taken to be 

f(t,x) = 2.04 sin(4707r/) 
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which models a periodic plane wave with a root mean square (rms) sound pressure level of 
117 dB. The frequency of 235 hertz is approximately halfway between the first and fourth 
natural frequencies of the system (as shown in [11], these occur at 65.9 hertz and 387.8 hertz, 
respectively). 

The dynamics of the uncontrolled system were approximated using 80 cavity basis functions 
(m x = m v = 8) and 11 beam basis functions (n = 12). The time interval of interest was taken 
to be [0, 16/235] which admitted 16 periods of the driving frequency, and time histories of the 
beam displacement at X = .3 and cavity pressure at X = .3,1^ = .1 on this temporal interval 
are plotted in Figure 7. 

The frequency plots of the uncontrolled beam displacement and cavity pressure in Figure 8 
exhibit not only the driving frequency but also transient responses at 65.9, 181.6,345.2,387.7 
and 519.5 hertz which are due to the natural frequencies of the coupled system (see [11] for a 
complete discussion of the dynamics and natural frequencies for the corresponding linearized 
system). In particular, the high energy response at 181.6 hertz indicates a strong excitation 
of the system at what corresponds to the frequency for the first mode of the uncoupled cavity 
(care must be taken when describing the dynamics of the system in terms of the undamped 
beam and cavity modes since the nonlinear coupling and beam damping yield system responses 
which differ somewhat from those of the isolated components). The presence of the multiple 
frequencies can also be seen in the time history plots of the uncontrolled beam displacement 
and cavity pressure in Figure 7. 

Control was then implemented by using Potter’s method to calculate the gains for the 
linearized system and feeding them back into the nonlinear system as discussed previously. 
The following results were obtained with an out-of-phase single pair (so as to create pure 
bending moments) of centered patches covering one half of the beam length as shown in Figure 
6. The quadratic cost functional parameters were taken to be d\ = d 2 = d 4 = 1,^3 = 10 4 
and R = 10~ 6 with d 3 chosen to have larger magnitude so as to more heavily penalize large* 
pressure variations. 

Figure 9 contains a plot of the controlling voltage u(t). As expected, it is periodic, and 
the magnitude remains below 25 V which is a physically reasonable voltage to apply to the 
piezoceramic patches. 

The application of the controlling voltage resulted in a high frequency transient response 
and 168 cavity basis functions ( m x = m y — 12) and 15 beam basis functions (n = 16) were 
needed to resolve the controlled system dynamics. 

From Figure 7, it can be seen that the controlled responses undergo a transient phase of 
approximately three periods and then are maintained at a low magnitude throughout the rest 
of the time interval. By calculating the rms pressure levels, it was determined that at the point 
(X, Y) = (.3, .1), the uncontrolled sound pressure level is 82.8 dB whereas the controlled sound 
pressure is reduced 15.7 dB to 67.1 dB. The level of reduction becomes even more significant 
as one moves deeper into the cavity since the strong cavity excitation in the uncontrolled case 
yields high magnitude pressure oscillations near the back wall which are uniformly reduced 
by the application of the controlling voltage. Finally, it is noted that the relative reduction 
in pressure is more significant than the reduction in beam displacement. This is due to the 
heavier penalization of pressure fluctuations through the choices cfo = 1 an d d 3 = 10 4 . 

The frequency plots of the controlled responses (in Figure 8) show that the dominant re- 
sponse is now at the driving frequency of 235 hertz. They also demonstrate the presence of 
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high frequency transient responses which are much more significant than those found in the 
uncontrolled case. This indicates that the interior pressure oscillations are reduced through 
two mechanisms when the controlling voltage is applied; the first is due to the reduced mag- 
nitude of the beam displacements while the second is due to the excitation of high frequency 
beam oscillations which couple less readily with the interior acoustic field. The combination 
of the two results in significantly reduced interior sound pressure levels. 




Fig. 7. Uncontrolled and controlled beam displacements and pressures at the points X = .3 and 
(-X\T) = (.3, .1) throughout the time interval [0, 16/235]. 
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5 Conclusion 


In this paper, we have discussed several of the issues which are involved in using piezoceramic 
patches as actuators in a nonlinear structural acoustics application. The patches affect the 
dynamics of the coupled system by contributing external forces and moments to the structure 
when a voltage is applied, and the first part of the discussion is centered around a description of 
the interactions between the patches and an Euler-Bernoulli beam and a thin cylindrical shell. 
In this discussion, care was taken to distinguish between the passive (material) contributions, 
due to the added thickness and differing material properties of the patch and bonding layer, 
and the active (external) contributions which result from the strains which are produced when 
a voltage is applied to the patches. 

As a result of the differing material properties and presence of the piezoceramic patches, 
the material and control parameters of the combined structure are piecewise constant in nature 
and hence lead to discontinuities in the moment and force resultants. This leads to difficulties 
in the strong form of the system equations when the moments are differentiated and is one 
motivation for using the weak or variational form where the derivatives are transferred onto 
the test functions. The weak form is also advantageous for many approximation schemes since 
it reduces the smoothness requirements for the basis elements. Finally, well-posedness issues 
were considered by posing the weak form in the context of sesquilinear forms. 

Due to the nonlinearities arising in the coupling between the beam vibrations and the 
interior acoustic field, LQR feedback control results could not be directly applied to the 
problem. Instead, gains corresponding to the linearized problem were calculated and fed back 
into the nonlinear system. As demonstrated by the results in the example as well as the more 
extensive set of examples in [14], this strategy is very effective for this problem. This is partly 
due to the weakness of the nonlinearity. By comparing the nonlinear results reported here 
and in [14] with the corresponding linear ones in [11], one can see that qualitatively, the two 
sets agree closely. This can be explained by the fact that the beam displacements are very 
small and hence the linearized coupling terms quite accurately approximate the true nonlinear 
expressions. 
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